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O ■ I- INTRODUCTION 

o ; 

' In recent years there has been much interest in braneworld scenarios 0, 0i 0] j where our four dimensional universe 

could be a brane embedded in a higher dimensional bulk spacetime. The basic ingredient of braneworld scenarios 

^ • is that while ordinary matter fields are confined to the brane, only gravitational interactions can propagate in the 

bulk. The second Randall-Sundrum model (RS II) ^ with one brane in an anti de Sitter (AdS) bulk is particularly 

Qv^ ' interesting in the point that, despite an infinite extra dimension, four dimensional general relativity (4D GR) can be 

recovered at low energies/long distances on the brane. Then, what is leading order corrections to the conventional 

T-H . gravitational theory? It was shown that the Newtonian potential in the RS II braneworld, including the correction 

K*" ■ due to the bulk gravitational effects with a precise numerical factor 0, j is given by 
I \ 

O^ '- , Gmm' ( 2f 

g: "(•■):=-— + 

' Here I is the curvature length of the AdS space and is experimentally constrained to be £ < 0.1 mm. It is natural to 
, ask next what are leading order corrections to cosmological perturbations. 

' Great efforts have beenpaid for the problem of calculating cosmological perturbations in braneworld scenarios (see, 
for example, Refs. 0, i, ll [IS [IE El El El El El El [S, 19, 20, 21, 22, 23, 24, 25, 26, [13). In order to correctly 
I • evaluate perturbations on the brane, we need to solve bulk perturbations, which reduces to a problem of solving partial 
differential equations with appropriate boundary conditions. Hence it is not as easy as in the standard four dimensional 
• • . cosmology. A pure de Sitter braneworld E3, and its variations such as "junction" models 0,0] and a special dilatonic 



. ^ ' braneworld [2J, |25|, |2y| are the rare examples where the bulk perturbations are analytically solved. In more generic 
, cases, one has to resort to a numerical calculation [TsL [2ll [2^. or some approximate methods [13. [lol [20l [23l| . There 
^ is another difficulty concerned with a physical, fundamental aspect of the problem; we do not know how to specify 
5^ , appropriate initial conditions for perturbations with bulk degrees of freedom. 

In this paper, we investigate cosmological tensor perturbations in the RS II braneworld, and try to evaluate leading 
order corrections to the 4D GR result analytically. For this purpose, we make use of the reduction scheme to a four 
dimensional effective equation which iteratively takes into account the effects of the bulk gravitational fields [2^ . As 
a first step, we concentrate on perturbations which is initially at a super-horizon scale. 

The paper is organized as follows. In Sec. II we define the setup of the problem discussed in the present paper 
more precisely. In Sec. HI we briefiy review the derivation of the effective equation of motion for tensor perturbations. 
In Sec. IV, using the long wavelength approximation we develop a method for evaluating leading order corrections, 
and evaluate the corrections considering slow-roll inflation and the radiation stage which follows after inflation. The 
shortcoming of the long wavelength approximation is cured by introducing another method suitable for analyzing the 
radiation phase in Sec. IV. Section V is devoted to summary and discussion. 
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II. SETUP OF THE PROBLEM AND BACKGROUND COSMOLOGICAL MODEL 



In this paper we investigate tensor perturbations in a Friedmann braneworld model with an warped (infinite) extra 
dimension 0, [s^, |^ IS E^l ■ To determine the spectrum of tensor perturbations observed at late time, we need 
to solve the perturbation equations with appropriate initial conditions. However, it is not an easy task to specify 
what initial conditions we should use in the context of braneworld cosmology. At an earlier stage of the evolution 
of the universe the physical wavelength of perturbation modes is much shorter than the bulk curvature scale. Then, 
the correction to the four dimensional standard evolution will not remain small. Therefore, in this regime it is not 
appropriate to discuss the dynamics of perturbation modes under the assumption that a modification to the four 
dimensional standard one is small. Hence our strategy is to isolate the issue of initial conditions from the rest. In 
this paper we discuss a relatively easy part, i.e., late time evolution after the effect of the fifth dimension becomes 
perturbative. 

Our discussion will be restricted to the perturbation modes whose wavelength is initially much longer than the 
Hubble horizon scale. We analytically derive a formula for small corrections to the solution for the tensor perturbation 
equation due to the presence of the fifth dimension. We mainly focus on the behavior of a growing solution, since 
a growing solution will be more important than a decaying one in general. Later we briefly mention the decaying 
solution in Sec. IVII 

The background spacetime that we consider is composed of a five dimensional AdS bulk, whose metric is given in 
Poincare coordinates as 

^2 

ds^ = — {dz^ - dt^ + Sijdx'-dx^) , 

with a Friedmann brane at z = z{t). Here £ is the curvature length of AdS space. In the original RS II model Q, 
a brane is placed at a position of fixed z, and Minkowski geometry is realized on the brane. The induced metric on 
z = z{t) is 

ds^ = a^it) - P)dt^ + Sijdx'dx^] , (1) 
a(t) = i/z{t), 



where the overdot denotes dt. From Eq. J^l, we see that the conformal time on the brane is given by drj = Vl — z'^dt. 
The brane motion is related to the energy density of matter localized on the brane by the modified Friedmann 
equation jHIMEllllIll as 



2a 



and 



i^---^, (3) 

where a = 3/47rG^^ is a tension of the brane. The quadratic term in p in the Friedmann equation modifies the 
standard cosmological expansion law. At low energies {£H <C 1), the brane motion is non-relativistic (i <C 1), and 
hence the conformal time on the brane rj almost coincides with the bulk conformal time t. 

As a background Friedmann brane model, we consider slow-roll inflation at low energies followed by a radiation 
dominant phase. During the slow-roll inflation the wavelength of perturbations stays outside the Hubble scale, and it 
re-enters the Hubble horizon only during the radiation era. 

We fix a model of slow-roll inflation by simply providing a function 

(4) 

where Ti. :— d,^a/ a = aH . During slow-roll inflation, e is supposed to be small. In the e limit, we recover the de 
Sitter case 

^ = 7, > Vl < Vo), (5) 



2??o - 7? ' 2770 - 77 ' 

where we have chosen the normalization of the scale factor and the origin of the ?7-coordinate so that a = aq and 
H = 1/770 at ry = r]o. 
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In the above alternatively we could specify the potential of the inflaton field localized on the brane. Then, for a given 
inflaton potential, the evolution of the background scale factor would differ from the conventional four dimensional one. 
However, as we will show below, the correction to the tensor perturbations coming from the slow-roll inflation phase 
is always small irrespective of the details of the evolution of the scale factor. Hence, we simply give time-dependence 
of the slow-roll parameter to specify a model. 

We assume that slow- roll inflation terminates at around rj = rjo, and a radiation era follows. During the radiation 
era the energy density decreases like p = pQ{ao/a)^ as usual. Then, solving Eq. |2Jl at low energies, we find that the 
scale factor behaves like 



(0) 



(6) 



with 



(0), , T] 



5a[ri) = aofH^ 



1 

24 



i (JL 

8 \vo 



(7) 



and Hq — {aorjo) ^. Correspondingly, we have 



(8) 



with 



(0) 

niv) 



6H{v) 



iv > m) 



(9) 



It is convenient to extend the definition of the slow-roll parameter e to the later epoch by 



e(?7) 



:=1-^. 



(10) 



III. TENSOR PERTURBATIONS ON A FRIEDMANN BRANE 



We use the method to reduce the five dimensional equation for tensor perturbations in a Friedmann braneworld at 
low energies to a four dimensional effective equation of motion, derived in Rcf. 28] . Here we summarize the derivation. 
Tensor perturbations on a Friedmann brane jSIiiliilil|l9l|20l|2ll|2^|23 are given by 



ds^ = — [dz^ - dt^ + (% + hij)dx'dx^] 



(11) 



We expand the perturbations by using ^^^(a;), a transverse traceless tensor harmonics with comoving wave number 
k . as hi 



— Yf'j{x)<^k{t^ z). Then the equation of motion for the tensor perturbations in the bulk is given by 

3 



-di 



-d. + d 



z 



(12) 



Hereafter we will discuss each Fourier mode separately, and we will abbreviate the subscript k. The general solution 
of this equation is 



dui^{uj)e-''^%pzfK2{pz) 



= / dLj'i'{LL;)e~ 



1 - 



{pzf (pzY 



16 



— 7 — In 



(f) 



(13) 



where p^ = —lo^ + fc^, and 7 is Euler's constant. We have chosen the branch cut of the modified Bessel function K2 
so that there is no incoming wave from past null infinity in the bulk. The coefficients ^(tj) are to be determined by 
the boundary condition on the brane fi''9^$|2=z(t) = 0, where n'' is an unit normal to the brane, or, equivalently, by 
the effective Einstein equations on the brane 34] , 



(4)G., = 87rGr„, -t-(8^G5 



E„ 



(14) 
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Here tt^^ is quadratic in the energy momentum tensor T^,y, and the first two terms on the right hand side are totally 
represented by the variables which reside on the brane. A projected Weyl tensor iJ^j^ :— Cafj_p^n"n^ represents 
the effects of the bulk gravitational fields, giving rise to corrections to four dimensional Einstein gravity in a fairly 
nontrivial way. In the present case, this can be written explicitly, and the effective four dimensional equation reduces 
to 



-2E 



(15) 
(16) 



where (j){t) := ^{t,z{t)) is the perturbation evaluated on the brane. At low energies {£H ^ 1), we can use the 
approximation like 



A:2)$, 



neglecting higher order corrections of 0{£*). Then, using the lower order equation to eliminate the higher derivative 
terms in —2E, we finally obtain our basic equation: 



(0) „ 



where 



2 (0)2 

5o[0] = 



£2 

1 



2nn') 



dujp (j)e 



In 



2a 



+ 7 



(17) 

(18) 
(19) 
(20) 



and the prime denotes dr^. The first term arises due to the non-standard cosmological expansion included in TL, while 
the last two terms are the corrections from the bulk effects E^i,. We can see that all terms are suppressed by £'^ (or 
£^ h\£) ^ . S2 is essentially nonlocal because of the presence of the log term. A time-domain expression for the quantity 
which appears in 52 can be rewritten, by eliminating the higher derivative terms in a similar way as before, as 



(21) 



where we have introduced a derivative operator 



IV. NEW LONG WAVELENGTH APPROXIMATION 



Let us begin with the case that the wavelength of perturbations is much longer than the Hubble scale. The long 
wavelength approximation not only allows us an easy treatment of perturbations, but also turns out to be sufficient 
to derive all the major corrections except for that coming from Sq. We would like to stress that what we call the long 
wavelength approximation here is different from a low energy expansion scheme or a gradient expansion scheme in 
the literature 



jproximatic 



General iteration scheme 



The equation we are considering takes the form of 



dl + 2nd^ 



-s\ 



(22) 



^ Hereafter, we will refer to both terms suppressed by £^ and by \nl as 0{fi). 
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This can be solved iteratively taking €^ as a small parameter. For this purpose, we write 



F( ~1 

'Piv) —4> iv)^ —4> iv) exp 



fiv')dv' 



(23) 



(0) 



A zeroth order solution ip by definition satisfies the equation obtained by setting the right hand side of Eq. H22|l to 
be zero. Then / obeys 



(0) 

dr) 4> 



(0) 



which can be integrated immediately to give 



^ (0) (0) 

dr]' (p S[(p ] 



(24) 



(25) 



Integrating this expression, we obtain the first order correction F. In the present case, S consists of three parts, 5o, 
Si and 52- We hereafter denote the corrections / and F coming from Si as fi and i^i, respectively. 



B. Long wavelength expansion of the zeroth order solution 



In the long wavelength approximation, we can write the zeroth order solution explicitly as 



drf I{ri') 



(0) 



~ -Akk'Iirj), 



where Ak is an amplitude of the perturbation which can depend on fc, and 



(26) 
(27) 

(28) 

^ — oo 

Here we simply keep the terms up to 0{k'^). 

We have already introduced a parameter e(ry) in Eq. H1U|) . which parametrizes the cosmological expansion. In the 
inflationary universe, it reduces to a slow-roll parameter and is assumed to be small. Slow-roll inflation assumes that 
time differentiation of e is 0{e^). Hence, we treat e as a constant during slow-roll inflation. The transition from 
inflation to the radiation stage occurs at around ij — rjo, and the behavior of e there depends on the details of the 
reheating process. During the radiation stage, e keeps again a constant value, e — 2, because the scale factor behaves 

like Woe rj. In general e reduces to a constant for the scale factor a'proportional to a power of the conformal time. 

(0) (0) 

Using e and the zeroth order solution in the long wavelength approximation, we can write iSi [</> ] and p as 



= H' 



9 (") 9 

Akk^n^ 



(0) (0) 

l-Hl-2eHl 



(29) 



2n' 



e' ~2e{l + e)n 



(0) 



2ek^ 



-2Akk^n^ { e'l -2e 



{l + e)nl-l |, (30) 

which are generic expressions not relying on any speciflc time dependence of the scale factor. 

Equation H30|l can be further rewritten into a more suggestive form. First, integrating the identity Ji) = 

(l + e)a2, it can be shown that 



(0) 



(0) 

■H_ 

(0)„ 



(0)2 , / 

e a drj . 
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Then, integration by parts leads 

(0) 

il + e)'-^!-!^^ r e'^S-Id^', (31) 

and thus we have 

e'l~2e^( e'^'cI^IdA. (32) 
a2 J -oo I 

Since e' is of second order in the slow-roll parameters during slow-roll inflation, we can safely neglect this term until 
the end of inflation. Moreover, for sufficiently smooth e this term is expected to be small, and e' = for the scale 
factor proportional to a power of the conformal time. Only at the time when the equation of state of the universe 
changes abruptly, i.e., at a sudden transition from inflation to the radiation stage, e' possibly becomes significantly 
large like a delta function. After such a violent transition, the integral in Eq. 1)32(1 leaves a constant value and thus 

this second term decreases in proportion to the damping factor / 



9 (0) o 



C. Slow-roll inflation 



Now let us consider slow-roll inflation on the brane. In this case Sq vanishes, since SH = by construction. In 

(0) , 

the previous section we have also shown that p (p vanishes at first order in the slow-roll parameters during inflation. 
Consequently, the only relevant term in the present situation is iSi . 
Equation together with Eq. ||2U tells us that 



Si 



and thus 



fiv) 



k'^e r , , c), 
d-q'en . 



(33) 



(34) 



Using the relations 7i ~a/aoJ7o and drj = ao?yo rf'a/ and neglecting the time variation of the slow-roll parameter, 
we obtain 



fiv) 



mo- 



ao ao 



(35) 



where = a{rii) and ry^ is the lower boundary of integration. Integrating this expression, we have 



fa,'' 



















(36) 



We chose the integration constant so that F vanishes in the limit of airj) — > oo. This means that we have renormalized 
Ak so that it represents the amplitude of fluctuations that we see at late times if slow-roll inflation lasted forever. 

Although our present long wavelength approximation is not valid when the wavelength Uik^^ becomes comparable 
or shorter than the Hubble scale H, we can still arrange to be sufficiently small so that we can neglect the first 
term in Eq. (|36|) . For such a^, we have an expression independent of a^, 



Fir,) 



2 1/2 



kH' 



(37) 



and we do not need to care about the choice of the "initial time" r/i . We see that the correction arising during inflation 

is very tiny, suppressed by the slow-roll parameter e in addition to the factor A:^^^/ a'^. For pure de Sitter inflation 
there is no correction, as is expected. 
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FIG. 1: Behavior of Si around the transition time with the vertical axis in an arbitrary unit. Red solid line shows the case 
of a sharp transition with s = 0.02ryo, while blue dashed line represents the case of a smooth transition with s — 0.5rjo. The 
wavelength of the mode is chosen to be krjo = 0.01. 
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FIG. 2: Behavior of S2 around the transition time with the vertical axis in an arbitrary unit. Red solid line shows the case of a 
sharp transition with s = 0.02r;o, while blue dashed line represents the case of a smooth transition with s = 0.6770. The choice 
of the parameter is the same as in Fig. Q krjo = 0.01. 

D. Transition to the radiation stage 

In this section, we investigate the effects of the transition from the inflation stage to the radiation stage. First, 
just for the illustrative purpose, we plot the behavior of Si and 52 for the modes well outside the horizon in the 
neighborhood of the transition time 770 in Fig. ^ and Fig|3 These plots arc for e given by 

e(r/) = tanh[(r/ - ryo)/s] + 1, (38) 

where the parameter s controls the smoothness of the transition. From the figures it can be seen that the corrections 
from E^^, Si and 52, become significant only around the transition time. 

From now on we shall consider the limiting case where e is given by e{'i]) ~ 29{ri ~ 779). Here we neglected the tiny 
effect of the non vanishing slow-roll parameter during infiation. For this instantaneous transition model, the scale 

__ . . (0) 

factor and the comoving Hubble parameter for rj < tjq are given by Eq. JSjl , while for rj > 7]q they are given by a and 
n in Eqs. iQ and ®. 
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Assuming such a step function-like transition, we can evaluate the corrections from <So at ?7 > rjo as 



9a2^o '^o 

kH^ .7] 



1 



1 + 2 



90al 



Vo 



(39) 
(40) 



Here the integration constant is chosen so that the correction vanishes before the transition. From this expression, 
we find that Fq diverges for a large r], and hence the correction to the final amplitude of fluctuations from So looks 
infinitely large. However, this is an artifact of the long wavelength approximation. At a late time during the radiation 
era the zeroth order solution becomes oscillatory. This oscillation suppresses the late time contribution. However, 
in the long wavelength approximation, this oscillation is not taken into account. In the succeeding section, we will 
develop another method which takes into account this oscillation. There we will find a finite correction. 
As for iSi, we have 



hiv) - - 



2pe r drf_ 

3 _ 9 
2 ~ 4 



1 + 5 



„/3 ; ' 



+ 



1 % 

•n 



V 



Long time after the transition (but of course before horizon re-entry), the correction becomes time independent. 



lim Fi{r]) 

77— >oo 



(41) 
(42) 

(43) 



The next step is to evaluate the correction coming from S2- This can be done as follows. First, for the instantaneous 



(0) 

transition <p becomes 



-^Akk^ 



ls{rj-m)-4e{v-vo)4 
m V . 

. (0) 



-4Akk^d^ 



Then, integrating by parts, the Fourier transform of (f) is obtained as 



P 



Vo 



TT J 7]' 



2Akk'' . 



(44) 



(45) 



with 

We separate <S2 into two parts, 



dujp 



(0) 



9(<^%) := 



dy- 



ciujrioiv-l) 



1 ( 



and, by substituting the expressions obtained above, we have 



Akk^ 
27r 



j du q{ur]o)dn 



+ 2Akk'^ 



7 + ] 



Inl* 



2ao 



, (0) 



V V 



where T := rj — Then, it can be integrated to give 

kH^ 



h{ri) = 



(0)„ 



+ 2 



1 
2^ 



1 , 

7 + In — 



duj qiojrio) 

M_ 
2ao 



ln(* 



In 



,2 \ 1 



(46) 



(47) 



.(T)^+/;^W-.o)^^}. (48) 
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The brunch cuts of log function and the path of w integration in the first term should be chosen so that the retarded 
boundary condition is ensured, which are presented in Fig. 13 The integration is dominated by « A:^, where we 
may approximate q{u}r]o) « q{kr]o) ~ 1/3 because krjQ is assumed to be small. Then, the first term can be written by 
using the formula presented in Appendix, and can be integrated to give a part of F2 (77) as 



3al 



cos[k{ri' ~ 770)] - 1 



+ d,, Mtj' - 770) ln(fc(77' - m))] + lS{7f - m) 



2ee 



2ee 

2,al 



2 / 2 

%ln(fc(77-77o)) + ln(fcr7o) 1-% 



[7 - 1 + ln(fc77o)] 



- 1 



In 



V - ??o 
V 



In (77/770 



+ 7 



where we have neglected the term proportional to cos[fc(7/' — 770)] — 1 in the integrand because this term brings the 
higher order contributions in k. The other part of i^2('7) is obtained by integrating the remaining parts in Eq. I|48l) as 



2k^e^ 



7 + In 



ke 

7 + In 



1 

5 
k£ 
2ao 



1 - 



1 

25 



4 
rf' 



- 1 



^ ln(77/??o) 



1 - 



1 

20 



1 - 



4 
25 



Combining the above two, we finally obtain the correction from 52, 



'f^OO flg 1^ 75 



2 

15 



7 + In 



/ k 



\aoHo 



-In 
5 



2 



(49) 



CORRECTIONS IN THE RADIATION STAGE 



So far we have worked in the long wavelength approximation, and found fc-dependent small corrections to the 
spectrum of tensor perturbations due to Si and 52, which correspond to E^j^i, term. It might be possible, however, 
that further corrections arise after the mode re-enters the horizon. We will show that such corrections are highly 
suppressed. A basic observation supporting this conclusion is that the contributions from term in the right 
hand side of Eq. H17|) become significantly large only at the transition time. It decreases in powers of 0(77) after the 
transition, and will be negligible when the long wavelength approximation brakes down. 

In the long wavelength approximation, we could not obtain a meaningful estimate for the leading order correction 
caused by unconventional cosmic expansion, Sq. Here we will give a more rigorous treatment, and will resolve this 
drawback of the approach taken in the preceding section. 

In the radiation stage. 



u{r]) := a{vi)cj){rii), 



(50) 



























1 ~^ 







FIG. 3: Branch cuts and the contour of the integration in Eq. 14811 for T > 0. For T < we close the contour on the upper 
half plane. 
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is a convenient variable. In terms of this new variable, we can rewrite the equation of motion (|17|l as 



u" + k^u : 



(51) 



with 



So = -^u. 



(52) 



If we neglect the deviation of the expansion law from the standard one, we have a oc 77 in the radiation stage and so 
a" — 0. This means that the first term on the right hand side gives a contribution of 0{£'^) from the modification of 
the expansion law as So in the preceding section. Thus, all the terms collected on the right hand side are of 0{£^), 
and at zeroth order u behaves just like a harmonic oscillator. The zeroth order solution is expressed by a linear 
combination of the standard modes: 



u {ri) = Ak 



-ikr] 



2k 



(53) 



where the coefficients a and f3 are determined by matching the solutions at the onset of the radiation stage. 
We define "energy" of the harmonic oscillator as 



fc2 



(54) 



Then, this energy is conserved at the zeroth order, i.e., 

5 =^(|ap + |/3n|AfeP= const. 

Taking into account the corrections of 0{P), the total variation of the energy can be estimated as 

2 2 



^ ^ /"OO 



dr], 



(55) 



with 



£' 



2a 
2a 



{u'ySo + c.c, 

{u'ySj + c.c, for j = 1,2. 



(56) 



The amplitude of the oscillator |v4i.|e'^[-'^l is proportional to a' /a. The factor a' /a here is to be attributed to 
unconventional cosmic expansion. Hence we have 



lim <^ 5R[Fo(7/)] 



77— >oo 



Ago 

(0) ' 

2£ 

lim 5R[F,-(^)] = 



for J = 1,2. 



(57) 



A. Corrections due to the unconventional expansion law 



We start with computing the correction due to the unconventional cosmic expansion, Fq, which we could not 
determine in the long wavelength approximation. To do so, we first need to determine the coefficients in Eq. H53|l by 
requiring that u*^"^ and its derivative are continuous across the transition from the de Sitter stage to the radiation 
stage. We assume for a moment that the perturbation in the initial phase is given by the de Sitter growing mode. 



v{r]) 



2?7o - V 



k\2r^o-vf 



(58) 
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Solving 



(0) , 



we obtain the coefficients as 



, (3^ a*. 



The energy of the harmonic osciUator at zeroth order is approximately given by 



(0) _ \Ak\^al 



Next we integrate 



(59) 



(60) 



(lap - - a(3*e-^''''> + a* Pe''"''>)\Ak\'' + c.c. 



The scale factor correct up to 0{£'^) is given in Eq. ©. Using this expression, we have 



a 
a 



2 776 



Then, integrating Eq. I|61(l . we have 



(0) 

2£ 



2a§ 

/2rr2 



a/3* / 



4fcr7o 
1 



^0 ^ 
2 



(61) 



(62) 



2 tr2 



1 ee 



- - ^^'yo - g(fc?7o) 



2i 



10 



(63) 



The first term in the last line is not suppressed for the fc — > limit. This is because this term arises due to the 
asymptotic behavior of da{ri). In fact one can see from Eq. 157|l that the expression for the physical amplitude of the 
perturbations does not have this term, and finally we have 



lim n[Fo{v)] 

77— ♦00 



1 k 



2 02 



10 al 



(64) 



B. Suppression at sub-horizon scales 



Here we shall show that corrections from the E^^.^ term are suppressed at sub-horizon scales and hence the result 
obtained in the preceding section by using the long wavelength approximation suffices for our purpose. 



Substituting =u j a into Si[(p], energy gain or loss due to this term is obtained as 



'-■1 



fc2 

„2 



/o), 



c.c. 



(65) 



Note that | 'tt'p and ')* it consist of just constant parts and oscillating parts with constant amplitudes. Therefore, 
it is manifest that when we integrate 5( the dominant contribution comes from around the lower boundary of the 
integral, 77 w 779, and so corrections arising in the sub- horizon regime are suppressed compared to those imprinted 
beforehand in the super-horizon regime. 
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A similar expression for £'2 can be obtained as 



8 a -I / a- 



2£2 

'2 



In 



kl 

(0, 

2 a: 



+ 7 



u'r + 



fc2 

„2 



+ C.C. 



(66) 



The second term is local and suppressed for large rj for the same reason above. To examine the late time behavior of 
the first term, we substitute the expression 153|l without limitation t] > rjo. This approximation is justified because 
the kernel J duje~^'^''^~^ ^p^ ln(p^/A:^) decays at least as fast as 1/(77 — 7/)^ for a large time separation, as will be shown 
in Appendix. Under this approximation for <f>{ri), it is easy to calculate its Fourier transform ^ 



(0) 



irio 1 



-{a + (3) for LO<-k, 



2ao ^/2fc 



X <^ (a ~ /3) 



for 



< fc. 



(a + /9) for LU > k. 



(67) 



The contour of the cj-integration in the first term of Eq. H66|l is such shown in Fig. and hence it can be separated 
into three parts as 



The first term in H66|l 



^ ^ 16 ^ ^ 



k—ioo 
k 



r-fe 



where 



ltk~ioo — i^- 



— k — ii 
-k 

i 



(a + /3). 



-ikrj 



ln(22/ - iy^) 



2k J- 



dyil-yYe'^'^^-yHl-y^). 



(68) 

(69) 
(70) 



We can see that all Is are suppressed for large values of krj, since their integrand becomes a product of a sooth 
function and a rapidly oscillating function. The dominant contribution in the integration of £' therefore comes from 
around the lower boundary of the integral. 

Based on the above analysis, we conclude that all the dominant corrections arise in the super-horizon regime, and 
thus it is sufficient to evaluate the corrections in the long wavelength approximation. 



VI. SUMMARY AND DISCUSSION 



We have investigated leading order corrections to the tensor perturbations in the RS II braneworld cosmology by 
using the perturbative expansion scheme of Ref. 28]. We have studied a model composed of slow-roll inflation on 
the brane, followed by a radiation dominant era. The unperturbed five dimensional bulk is AdS space with curvature 
radius £. In our expansion scheme the asymptotic boundary conditions in the bulk are imposed by choosing outgoing 
solutions of bulk perturbations, whose general expression is known in the Poincare coordinate system. Hence, the 
issue of bulk boundary conditions is handled without introducing an artificial regulator brane. This is one of the 
notable advantages of the present scheme. 

We set the initial condition when the wave length ak~^ is already longer than the Hubble radius during slow-roll 
inflation. We will not lose much by neglecting the modes which are already inside the Hubble scale at the time of 
transition to the radiation era, since they are not cosmologically so interesting. We compared the resulting amplitudes 
of fluctuations after the horizon re-entry with and without the effect of an extra dimension. To do so, we normalized 
the amplitude so that the late time amplitude of fluctuations becomes identical in two cases if the slow-roll inflation 
lasts forever. As the reference without the effect of an extra dimension, we took a growing mode solution. The 
corrections due to gravity propagation through the fifth dimension dominantly come from the contribution around 
the transition time. Therefore they can be estimated by using the long wavelength approximation with a sufficient 
accuracy. Since the correction due to the modified cosmic expansion comes from relatively late epoch, it is necessary 
to take into account the oscillatory behavior of the solution after the re-entry to the Hubble horizon. 



^ Depending on the choice of the integration path near rj = 0, the Fourier transform of 4>{r^) changes by a constant independent of uj. Here 
we took the principal values. 
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Combining the results obtained in Eqs. H43() . (|49|l and (|64() . we find that the amplitude of a fluctuation with 
comoving wave number k is modified by a factor e^^^' due to the effect of an extra dimension with 



lim ^[F{f])] 

77— >Ctt 



13 
150 



4 
15 



7 + In 



/ k 



\aoHo 



5 V 2 



(71) 



where oq and Hq are the scale factor and the Hubble parameter at the transition time, and 7 is the Euler's con- 
stant. Leading corrections are proportional to k'^P /a^ or k^i^ / a^\ogk, as is expected from the dimensional analysis. 
However, our calculation here determined the precise numerical factors analytically. 

Throughout this paper, we have dropped the contribution from the decaying mode for simplicity. Since the decaying 
mode during the initial slow-roll infiation phase is usually suppressed when the wavelength is longer than the Hubble 
scale, it does not give a significant effect. However, in the context of braneworld, we have not yet understood how 
to fix the initial conditions for fluctuations. Therefore, in principle, there is a possibility that the contamination of 
the decaying mode can be extremely large, although it seems quite unlikely. In such a case the contribution from the 
decaying component cannot be neglected. 

It is quite easy to estimate the leading order correction due to the decaying mode, since it does not pick up the 
effect of the presence of an extra dimension. Here the effect is not due to the non-trivial evolution of a solution but 
totally due to unconventional initial condition. Adding a decaying mode will modify v given in Eq. ()58|l to 



v{ri) 



2?7o - V 



1 



fc2(2ryo-ry)2 



where Cd is a complex number which parametrizes the amplitude of the decaying component. Repeating a similar 
calculation leading to Eq. H64|l . we find that the relative change in the amplitude due to the decaying mode is given 

by 



"fiiCdV-P^ikvor. 



(72) 



Thus the correction due to the decaying component is more suppressed in the sense of the power of k, but it might 
be more important than the other corrections for fcryo > 5R[Cd]~^ if 5i[Cd] is extremely large. 

Lastly we would like to mention the limitation of the present work and some future issues. The perturbative 
expansion scheme that we have adopted in this paper is valid for scales larger than the AdS curvature length {k£/a <ti 1) 
at low energies {£H <^ 1); besides our study is restricted to initially super-horizon perturbations. We gave an initial 
condition for a solution by hand, but of course such an approach is not satisfactory. To determine the initial condition, 
we need to investigate the evolution of perturbations at the initial phase where the wavelength of the mode is too 
short and reduction to an effective four dimensional problem is not possible. Furthermore, to seek for interesting 
effects that might arise at high energies {£H > 1), we have to develop a new formalism, which is a next challenge and 
will be hopefully reported in our future publication. 
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APPENDIX A: DETAILS OF CALCULATIONS 



In this appendix, we will show the following formula: 



duj' 



-iujT 1 



cos(A:T) - 1 



- dr [0(T) In(fcr)] + 7<5(r) 



(Al) 



where — k^ — uj'^ and 7 is Euler's constant. The branch cuts are all running on the lower half complex plane of uj 
(FigO, which ensures the retarded boundary conditions. 
We divide the above integration into two parts: 



duje 



In 



duje-"^'^ In 



(A2) 
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The first term of the right hand side can be rewritten as 

J duje-"^'^\n (-^^ = J dioe-"^^ [ln(w - k) + \n{uj + k) - 21nw 
The integration can be performed as 

j dwe-"^^ - k) = 9{T) f {~i)dre-^'^'^^-''-^ \n{ 



) + / {-i)dre~ 





T(k- 



ln(r 



~2Tr9{T)- 



ikT 



T ' 

and the other two integrals are obtained by setting k 



diue-'^^^-^l ^ 



k and fc — > 0. Thus we have 
cos(fcr) - 1 



T 



This expression is regular at T = 0. To evaluate the second term in the right hand side of Eq. ljA2|l . we need a trick 
because the integrand suffers from a slow fall-off at large values of It should be interpreted as 



2idj 



-iujT 



duj- 



■In 



2idT { e{T) 



dX— (In(-A) - ln(/jT)) 
A 



(A3) 



where we put A = iluT. Now the branch cut and the contour of the integration on A plane are such shown in Fig2| 
Thanks to this trick, the expression for the second term becomes a well-defined one, though the result does not change 
as long as convergence is guaranteed. Then, performing the integration as 



J (In(-A) - In(fcr)) = liind^ J dXe-^i-X)"-^ + 2TTiln{kT) 

= lim [(e-""^ - e''"')T{i^)] + 2m ln{kT) 
= 2TTi [In(fcT) + 7] , 



we finally obtain the formula (|A1(I . 

For T > the formula ljAl|l reduces to / dcoe"^'^'^ In (p^/fc^) — —4tt cos{kT)/T and hence the kernel that appears 
in ^2 becomes 



due 



327r 



y3 



3 \ 

\os{kT)-—^s\n{kT) 



for T > 0. 



(A4) 



From this expression we find that the non-local source 52 does not keep the past history for a long time. 





.. > 




— ^ 




■ 4 



FIG. 4: Branch cut and the contour of the integration in Eq. IIA3l l. 
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